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We study the finite-temperature phase structure and the transition temperature of QCD with 
two flavors of dynamical quarks on a lattice with the temporal size N t = 4, using a renormalization 
group improved gauge action and the Wilson quark action improved by the clover term. The region 
of a parity-broken phase is identified, and the finite-temperature transition line is located on a two- 
dimensional parameter space of the coupling (/3 = Q/g 2 ) and hopping parameter K. Near the chiral 
transition point, defined as the crossing point of the critical line of the vanishing pion mass and 
the line of finite-temperature transition, the system exhibits behavior well described by the scaling 
exponents of the three-dimensional 0(4) spin model. This indicates a second-order chiral transition 
in the continuum limit. The transition temperature in the chiral limit is estimated to be T c = 171(4) 
MeV. 
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I. INTRODUCTION 



With the development of recent experimental projects toward a detection of the high-temperature quark-gluon 
plasma in relativistic heavy ion collisions, it is urgent to theoretically establish thermal properties of QCD directly 
from its first principles. 

For the pure gluon system, numerical simulations on the lattice are quite advanced by now, with the critical 
temperature and the equation of state, which are the basic thermodynamic quantities needed for phenomenological 
studies, well established in the continuum limit [EJ. Full lattice QCD simulations incorporating dynamical quarks, 
however, have not reached such a status. In particular studies on spatially large lattices, which are required for 
determination of equation of state, are still limited to lattices of a small temporal size, N t = 4 and 6 0], due to large 
computational power needed for dynamical quark simulations. 

A possible way to bring the simulation close to the continuum is to employ improved lattice actions: the continuum 
limit may be taken for thermodynamic quantities from coarse lattice spacings and correspondingly small temporal 
lattice sizes such as Nt = 4 and 6 or even smaller. For the pure gluon system, there are a number of tests of this idea 
reported some time ago || . More recently it is shown Q that the equation of state calculated for a renormalization- 
group (RG) improved gluon action || agrees with that for the plaquette action || in the continuum limit. In the 
present work we pursue an extension of this pure gluon study to full QCD with two flavors of dynamical quarks [j. 

Here we focus our attention to the Wilson-type quark action, and study the effect of improvement for both gluon 
and quark sectors. Compared to Kogut-Susskind (staggered) quark action ||, a clear advantage of the Wilson- type 
quark actions is manifest flavor symmetry. On the other hand, chiral symmetry is explicitly broken causing some 
subtleties in the analysis. This leads to large lattice artifacts at the range of lattice spacing a -1 « 1-2 GeV extensively 
used in finite temperature simulations, when the combination of the plaquette gluon action and the standard Wilson 
quark action is used. The most notable artifact is an unexpected strengthening of the finite-temperature transition at 
intermediate values of the quark mass |9|Jl0(| . This undesired behavior disappears when an RG- improved gauge action 
is employed pT| . A comparison with other action combinations, the plaquette gluon action and a clover- imp roved |l2j ] 
Wilson quark action p3|, or the Symanzik-improved gluon action and a clover- improved quark action flTif , indicates 
that improvement of the gluon action is essential in removing the dominant lattice artifacts at finite lattice spacings. 
On the other hand, we expect that improvement for the quark action renders the behavior of fermionic thermodynamic 
quantities close to continuum. 

Previous studies with the standard Wilson action have shown that the phase diagram in the plane of coupling and 
hopping parameter (fl = 6/g 2 : K) contains a phase with spontaneous breakdown of parity- flavor symmetry |p^| . The 
boundary of this phase is characterized by vanishing of pion (screening) mass. The interplay of the boundary with 
the line of finite-temperature transition plays a crucial role in understanding the finite-temperature transition with 
Wilson-type quark actions fll6| , [To 17|. Thus, prior to calculation of thermodynamics quantities, an examination of the 
finite-temperature phase structure is indispensable when one uses the clover-improved Wilson quark action coupled 
to the RG-improved gluon action. 

A particularly important point concerns the universality property of the finite-temperature transition. For the 
standard staggered quark action, extensive scaling analyses carried out for Nt = 4 failed to establish either the 
0(4) scaling expected from the effective sigma model analysis or the 0(2) scaling suggested from the actual 
symmetry of the staggered quark action p9|-p2]|. On the other hand, a study with the standard Wilson quark action 
coupled to the RG-improved gauge action found that the 0(4) scaling ansatz describes well the chiral condensate data 
^T| , P7| . Our interest here is to explore whether the 0(4) scaling is maintained for the clover- improved Wilson action. 
Such an analysis is also necessary to quantitatively locate the point of transition for massless quark. Combining 
this information with measurements of hadron masses at zero temperature, we obtain an estimate of the critical 
temperature. 

In this paper we adopt specifically the improved gluon and quark actions which is resulted from a systematic 
comparative study of various combinations of improved actions ]25| ] . The combination is the same as that we have 
used in our zero-temperature studies reported elsewhere [ pi|j2l| . 

This paper is organized as follows. The action and the simulation parameters are described in Sec. ||. The phase 
structure for our action is studied in Sec. [II. Section IV is devoted to the issue of the 0(4) universality. Assuming 
the 0(4) critical exponents, we then extrapolate the transition line to obtain a precise estimate of the chiral tra nsiti on 
point. The transition temperature in physical units is discussed in Sec. VI. A brief conclusion is given in Sec. VII. 



II. SIMULATION 
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A. Choice of the action 



We employ the RG-improved gluon action || 

S g = -/?{ co ]T W^\x) + Cl £ W^(x)}, (1) 

with /3 = 6/# 2 and ci = —0.331 and c = 1 — 8ci, coupled with the clover-improved Wilson quark action defined 
by 

• V '« y ^ T L l) -;'l, r (2) 

fx,!/ = fcy - - 7 M )£/x,ji&c+£,3/ + (1 + 7f*)E4 lA1 $c,v+/i} _ S xy CswK ^ o^F^. (3) 

Here i 7 ^,, is the lattice discretization of the field strength, 

F^v ^t(/^ii/ ~ (4) 
with f^ v the standard clover-shaped combination of gauge links, and we adopt the meanfield-improved clover coefficient 

c sw = (W lxl )~ 3/i = (1 - 0.8412/T 1 )- 3 / 4 , (5) 

using VF lxl calculated in the one-loop perturbation theory l|. This choice of Csy/i when expanded in /3, agrees well 
with the actual one-loop result csw — 1 + 0.678(18)//?+- • • pq| . Furthermore, the one-loop values of W^ lxl reproduce 
the mean values of plaquette within 8% for the range of (3 and K of our studies [^,^5) . 

This combination of improved gauge and quark actions was tested in our full QCD comparative study [§3|. We found 
that the scaling violation in light hadron masses and the violation of rotational invariance in a static quark potential 
are both small with this combination of actions already at a -1 w 1 GeV where our finite temperature simulations are 
made, as compared to a _1 > 2 GeV needed for the standard plaquette gauge and Wilson quark actions. 

B. Details of simulations 

We make our finite-temperature simulations on a lattice with a temporal extension N t = 4. Zero-temperature 
simulations with temporal sizes comparable or larger than spatial sizes are also made to calculate hadron masses for 
the purpose of fixing the scale. 

Our simulation parameters are summarized in Tables [ifV. Runs are carried out on lattices of size 16 3 x 4 and 



16 4 at 13 = 1.8-2.2, and of size 8 J x4 and 12 3 x 24 at j3 = 1.2-1.7. Except for the case of the 12 3 x 24 lattice, we 
use anti-periodic boundary conditions in the temporal direction for quarks. All other boundary conditions are taken 
to be periodic. Simulations are carried out at seven to eleven values of the hopping parameter K for each 0, which 
correspond to mps/my ~ 0.4-0.9 for the 8 3 x 4 lattices, and 0.6-1.0 for the 16 3 x 4 lattices. 

The Hybrid Monte Carlo algorithm is employed to generate full QCD configurations with two flavors of dynamical 
quarks. Details of the simulation are basically the same as in our zero temperature studies pi] , p5|] . The molecular 
dynamics time step St is chosen to yield an acceptance rate greater than about 80%. The length of one trajectory 
is unity in most cases, but is reduced down to 0.25 at several points close to the critical line to keep the acceptance 
rate of about 80% within 200 molecular dynamics time steps. The inversion of the quark matrix is made with the 
BiConjugate Gradient Stabilized (BiCGStab) method for K < K c , where K c is the critical hopping parameter at 
zero temperature where pion mass vanishes. For K > K c , we use the Conjugate Gradient (CG) method since the 
clover- Wilson operator may have negative eigenvalues, in which case the BiCGStab algorithm fails to converge. 

We measure gluonic observables, Wilson loops and the Polyakov line, at every trajectory. Hadron (screening) 
propagators are calculated at every 2-5 trajectories using both a point source and an exponentially smeared source 
for quark, and a point sink. For hadronic measurements on the 8 3 x 4 lattice, we double periodically the lattice in 
one of the spatial directions. On the 8 3 x 4 and 12 3 x 24 lattices, we determine hadron masses by a single hyperbolic 
cosine fit to the propagator calculated with the smeared source, because of its clear plateaus of effective mass. On 
the 16 3 x 4 and 16 4 lattices, since a plateau is sometimes less clear with our smeared source, we perform a combined 
fit using both point and smeared sources. 
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We measure the current quark mass denned through an axial vector Ward-Takahashi identity p7| , pq] , 

V M A M = 2m q P + 0(a) (6) 

where P is the pseudoscalar density and the //-th component of the local axial vector current. In practice we 
make a simultaneous fit of the two-point functions (A fJi (t)P(Q)) and (P(i)P(O)) to extract the pion mass m, and the 
amplitudes (0\P\ir(p = 0)), (Q\A^\ir(p = 0)}, from which we compute 

We choose \i = 4 for zero temperature simulations, and fj, = 3 at finite temperatures for which the screening masses 
are determined along the z axis. We also use an alternative definition of quark mass given by 

m q = (l/2a) (K- 1 - K- 1 ). (8) 

While these different definitions of m q give different values at finite /3, they converge to the same value in the continuum 
limit p^ , p5[ | . In this work, we concentrate on the phase structure and critical properties around the finite temperature 
chiral transition point. Therefore, in this work, we ignore the renormalization factors Zp, Za and Z m since they are 
regular around the transition point. 

We use the zero-temperature vector meson mass my as a measure of the lattice scale. For (3 > 1.8, detailed hadron 
mass data for light quarks are available from our extensive calculations |^4|. We use the results of hadron masses 
obtained in this reference for our analysis. These data, however, do not fully cover the region of heavy quark where 
the finite-temperature transition is located for Nt = 4 at (3 > 1.8. We evaluate hadron masses in these regions, and 
interpolate them by a cubic spline formula to the point of transition. 

No previous data are available in the region of strong coupling (3 = 1.3-1.7. Thus the hadron masses in this region 
are measured in the present work. We fit the vector meson mass in terms of pseudoscalar meson mass mps by an 
ansatz inspired by chiral perturbation theory, 

niya = Ay + By(mpsa) 2 + Cv(mpsa) 3 ■ (9) 

The fit curves are illustrated in Fig. |l|. The lattice spacing defined by identifying my(T = 0) = m p — 770 MeV in the 
limit of zero pion mass is given in Table 0. 

Errors are determined by a jack-knife method. From a study of bin-size dependence, we adopt the bin size of 1-10 
configurations (i.e. 2-50 trajectories) for hadron masses and 5-20 trajectories for other gluonic quantities. 



III. PHASE STRUCTURE 

The phase diagram obtained in our study is summarized in Fig. ||. Details are discussed in the following subsections. 
Thin open symbols show the simulation points on the Nt — 4 lattices. The solid line connecting symbols denoted 
as il K c (T — 0)" represents the critical line of vanishing pion mass at zero temperature. Other objects on the phase 
diagram are obtained on the Nt = 4 lattices. The K t -\me shows the location of the finite-temperature transition, 
above which the system is in the high-temperature phase. We expect the i^-line to cross the K C (T = 0) line. The 
crossing point is a natural candidate to be identified as the point of chiral transition jl0| . The shaded region shows 
the parity-broken phase p5[ . The lower boundary of the parity-broken phase lies closely above the K C (T — 0) line. 



A. Critical line at zero temperature 

The location of the critical line at zero temperature K C (T = 0) where the pseudoscalar meson mass vanishes is 
essential information for discussion of chiral properties with Wilson-type quark actions. To determine K C (T — 0), we 
extrapolate rrip S as a function of l/K using a quadratic ansatz: 

<( m p Sa ) 2 ^Ps(l-^) + Cps(l-i-) 2 . (10) 

The results for mp s are illustrated in Fig. || in the strong-coupling region of /3 = 1.3 — 1.7. Our values for K C (T = 0) 
are summarized in Table [v] and plotted in Fig. |2[ 
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Towards the weak coupling limit, = oo, the K C (T = 0)-line gradually approaches the free Wilson quark value 
1/8. For the case of the standard unimproved Wilson quark action, K C (T = 0) is a monotonically decreasing curve 
connecting K w 1/4 at = and 1/8 at = oo. For the case of our clover-improved quark action, the K C (T = 0)-line 
shows a maximum in K at « 1.47, and decreases as we lower /? below this value. This is likely to arise from our 
choice (p|) for the clover coefficient, which diverges as approaches 0.8412. 



B. Finite-temperature transition 



We identify the finite temperature transition point K t from inspection of the Polyakov line and Wilson loops. In 
Figs. H and|| the Polyakov line and its susceptibility at = 1.8-2.2 on the 16 3 x 4 lattice are shown. We fit the peak 
of the susceptibility by a Gaussian form using 3 or 4 points near the peak, except for = 2.2. The results for K t from 
the fit are given in Table [v|. The ratios mps/my and T pc /mv with zero-temperature meson masses interpolated to 
the K t point are also summarized in the table, where T pc is the pseudocritical temperature at K t . These results are 
used in Sees. and [v| 

On the K t -\ine thus determined, other physical observables, such as the plaquette, also show rapid changes. Smooth- 
ness of the data around K t suggests, however, that, for the range of quark mass we studied, the finite-temperature 
transition at K t is an analytic crossover. From Fig. |^ and similar plots for other observables, we see that the crossover 
becomes monotonically weaker with increasing 0. Increasing on the K t line corresponds to increasing the distance 
to the critical line, i.e. increasing the quark mass. In the limit of infinite quark mass K — 0, the K t line will end at 
the first-order deconfinement transition of pure gauge theory, which is located at = 2.2863(10) for N t — 4 |]30|fl. 

For the smaller region [0 < 1.7) where simulations are made on an 8 3 x 4 lattice, data for the Polyakov line 
become too noisy to determine the position of K t . We alternatively take the spatial plaquette shown in Fig. |, and 
identify K t from its rapid changes. The range of K t for = 1.6-1.7 thus estimated is summarized in Table VI. 

The JsT t -line approaches the zero-temperature critical line K C (T — 0) as decreases (see Fig. ^). We expect the two 
lines to cross at a point (0 c t, K c tj- This point is a natural candidate for the point of chiral transition since the pion 
mass, and hence also the quark mass defined through the Ward-Takahashi identity, vanishes at zero temperature Jl0| . 
The location of this point obtained by scaling analyses (see Sec. |v|) is marked by a star in Fig. [| 

Our estimate for K t stops at = 1.6. For 0<1,5, because the parity-broken phase appears at intermediate K as 
discussed in the next subsection, a more detailed study is required to identify the location of K t . 



C. Parity-broken phase 



For Wilson- type quark actions, chiral symmetry is explicitly broken by the Wilson term away from the continuum 
limit. The appearance of massless pseudoscalar meson (pion) at K C {T — 0) is, therefore, not due to spontaneous 
breakdown of chiral symmetry. Analytic calculations in the strong coupling limit and numerical simulations at 
intermediate couplings show that there exists a region in the phase diagram in which parity-flavor symmetry is 
spontaneously broken ("parity-broken phase") |l5| , [lq| . In this picture, the pion is understood as the zero mode 
of a second-order phase transition that takes place along the K C (T = 0) line, signaling spontaneous breakdown of 
parity-flavor symmetry. 

At zero temperature, the parity-broken phase extends from the strong coupling region towards the weak coupling 
region, forming sharp cusps which touch the weak coupling limit = oo at K = 1/8 and at four other values of K. 
In the present work, we concentrate on the branch of the parity-broken phase with the smallest positive values of K. 
The lower boundary of this branch is the usual critical line denoted by K C (T = 0) in Fig. ||. 

For finite temporal size Nt, the cusp of the parity-broken phase retracts from the weak coupling limit to a finite 
value of 0. Let us denote the position of the cusp as (/3 C usp, -Kcusp), and the boundary of the parity-broken phase as 
K c {N t ). This boundary is double- valued for < cusp , whose two branches we denote as K l ° wel (N t ) < K^ pcl (N t ). 

In the low temperature phase K < K t , we expect a massless pion to appear as K is increased. Therefore, the lower 
part of the if c (AT t )-lirie should be located near the zero-temperature critical line K C (T = 0), with the difference being 
0(a). On the other hand, we do not expect a massless pion in the high temperature phase K > K t . Hence the whole 
of K c (N t )-\me should be in the low temperature phase. These considerations imply that the i^ t -line runs closely past 
the cusp of the K c {N t ) line. Therefore we expect that the the chiral transition point is located near and above the 
cusp M: 

0ct=0cu Sp + O(a), (11) 
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with /3 c t _ /3 C usp- The possibility that the chiral transition point actually agrees with the cusp is not excluded. For 
the case of the standard Wilson quark action coupled with the plaquette gauge action, (3 c t has indeed been found to 
be close to /3 C usp [^6|- Similar results were obtained for the case of the RG- improved gauge action with the standard 
Wilson quark action [17] . 

Figure ^ shows results for pion screening mass squared (mpga) 2 with our action on an Nt = 4 lattice. Solid lines 
are linear or quadratic fits in 1/K. In the upper figure showing results for 1.7 > j3 > 1.6, the two fit lines for (3 = 1.6, 
from large or small values of 1/K, cross each other. This indicates that mps remains finite for all values of K, and 
hence the parity-broken phase does not exist for (3 > 1.6. 

In contrast, when we decrease (3 down to 1.5 and 1.4, we see evidence for two values of K c (N t = 4) (lower figure in 
Fig. [t]) corresponding to the upper and lower boundaries of the cusp of the parity-broken phase. 

We determine the location of K c (N t — 4) by an extrapolation of (mpsa) 2 linearly in 1/K, using the lightest two 
to three points. We confirm that a quadratic extrapolation gives consistent results. Results of K c (N t — 4) are 
summarized in Table [y|. We find the gap between two values of l/K c (N t ~ 4) to be 

A(l/K c ) = 0.020(12) at £=1.5 (12) 
= 0.159(12) at £=1.4. (13) 

For the upper K c (Nt — 4) line, we also try to estimate its location by extrapolating (mpsa) 2 linearly in /3, at fixed 
K = 0.16, 0.1625, 0.165, and 0.17, as shown in Fig. @. Results are listed in Table |v|. As is seen in Fig. ||, the points 
for K c (N t = 4), obtained from the 1/K fits and the (3 fits, together form a smooth curve. 

An enlargement of the phase diagram around the parity-broken phase is given in Fig. ^. Since the parity-broken 
phase is absent at (3 — 1.6 the boundaries of the parity-broken phase should terminate at a cusp at (3 — 1.5-1.6. 
In order to determine the location of the cusp point more precisely, we fit the upper and lower K c {N t = 4) data 
separately by a quadratic ansatz. From these fits, we obtain 

/3 cusp = 1.538(46), ^ cusp = 0.15687(69), (14) 

where the errors are estimated from the crossing points of the one standard deviation error bands of the two fits. 
From the discussions around (^Tj), this value for /3 cusp gives a lower bound of (3 ct . 



IV. 0(4) SCALING OF CHIRAL CONDENSATE 



From universality arguments, the finite-temperature QCD transition near the chiral limit is expected to be described 
by an effective a model fT§| ]. According to this description, the transition for two quark flavors in the chiral limit 
is either first order or second order depending on the strength of anomalous coupling which breaks axial Ua(1) 
symmetry. When this effect is negligible, the transition is of second order, and one expects that the critical properties 
are described by those of the three-dimensional 0(4) Heisenberg model. Hence tests of 0(4) scaling provides us with 
a useful way to study universality properties of the chiral transition for two-flavor QCD Jl9[ |. 

In the 0(4) Heisenberg model, the order parameter is given by the magnetization M. Near the second order 
transition point, M satisfies the following scaling relation; 

M/h 1/s = f{t/h 1/f3S ). (15) 

where h is the external magnetic field, t = [T — T c (h = 0)]/T c (h = 0) is the reduced temperature, the exponents have 
the values UBS = 0.537(7) and 1/5 = 0.2061(9) Q, and the 0(4) scaling function f(x) is also known Q. 

In Ref. JUL it was shown that the 0(4) scaling of ( |l5|) is well satisfied for the standard Wilson quark action 
combined with the RG-improved gauge action under the identification t ~ (3 — (3 c t, h ~ m q a, and M ~ (ipip), where 
quark mass defined by the axial vector Ward-Takahashi identity was employed. It is important to note that the naive 
definition of the chiral condensate {ipi^} is not adequate for Wilson-type quarks because chiral symmetry is explicitly 
broken. A proper subtraction and renormalization are required. A properly subtracted (ipip) can be defined via an 
axial Ward-Takahashi identity |^8| as first employed in Ref. |llj : 

(^) S ub = 2m q aZ^2(Tr(x)n(0)) (16) 

X 

where ■n(x) = q x "/c,q x with q x the lattice quark field. For the normalization coefficient Z we adopt the tree value 
Z = (2K) 2 , which is sufficient for our study of critical properties because Z is regular at the finite temperature 
transition point. 
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Figure [h] shows the results of ('ipip)sub as a function of Ward-Takahashi identity quark mass from our action at 
P = 1.8-2.2 obtained on a 16 3 x 4 lattice. We perform a fit to the 0(4) scaling function 0j by adjusting the value 
of j3 c t as well as the scales for t and h, with the exponents fixed to the 0(4) values Since scaling is expected 

toward the chiral transition point corresponding to massless quark, we start our fit with the entire data set, gradually 
reducing the range of (P,m q a) toward smaller (3 and m q a until an acceptable x 2 /Ndf is obtained. We find this 
condition to become satisfied for the range (3 — 1.8-1.95 and 2m q a = 0.0-0.9 with x 2 /Nof — 0.82 for Nof — 33, for 
which we obtain 



p ct = 1.469(73). (17) 

This fit is shown in Fig. fy] 

In this test we use results for m q a obtained on the finite-temperature 16 3 x 4 lattice. We repeat the test using m q a 
determined on the zero temperature lattice, 16 , at the same values of (P,K). Results for quark mass extracted on a 
16 4 and a 16 3 x 4 lattice are mutually consistent for the range (3 = 1.8-1.95 and 2m q a = 0.0-0.9. As shown in Fig. [IJ, 
we again find a good agreement with the 0(4) scaling function with x 2 /Ndf = 1-18 for Nof = 22 at 

p ct = 1.462(66). (18) 

The number of data points are smaller due to fewer simulation points on the 16 4 lattice. 

The good consistency of our condensate data with the 0(4) scaling suggests that the chiral phase transition of 
two-flavor QCD in the continuum limit is of second order. The estimated value for (3 c t is somewhat low compared to 
the result for the cusp (3 cusp = 1.538(46) obtained in Sec. |III C| . The reason for this trend is not clear at present. A 



possible origin is lack of condensate data below (3 = 1.8, closer toward the chiral transition point, in the scaling fit. 
We leave an examination of this point for future work. 



V. SCALING OF THE PSEUDO-CRITICAL TRANSITION POINT 



Let us denote by t pc the pseudo-critical point defined as the peak position of the magnetic susceptibility at finite 

h. From the scaling relation (|l^), we expect 

t pc cx h 1 " 35 (19) 
to be satisfied for this quantity. Identifying t pc = (3 pc — Pct, where (3 pc is the /3-coordinate of the i^-line, we expect 

P P c = Pa + Bph 1 ^ 6 (1 + 0(h)) (20) 

near the chiral transition point. 

In Fig. |l3] we plot the pseudo critical point P pc choosing 

h = m q a (21) 

i. e., the quark mass defined by the axial vector Ward-Takahashi identity. Errors attached to m q a include the 
uncertainty in the location of the K t line for each value of P pc , which actually dominates the error. The precision 
of our results is not sufficient to attempt a fit taking the exponent 1/PS as a free parameter. We therefore make a 
fit adopting the 0(4) value for this exponent. Taking 1.6 < P pc < 1.95 for the fit range, since this was the range 
acceptable for the scaling analysis of the chiral condensate, we obtain 

p ct = 1.557(28) (22) 

with x 2 /N D f = 1-04 for N DF = 6. 

In Fig. ^ we show results of a similar analysis, making an alternative choice, 

h = 1/K t - 1/K C , (23) 

for the external field. The fit using this identification leads to 

p ct = 1.613(17) (24) 

for data in the same range 1.6 < P < 1.95 with x 2 I^df = 1-60 for Nof = 6. 
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As the third possibility, we consider 

h = (mps/my) 2 (25) 

obtained on the zero-temperature lattice. This choice has the nice feature that the location of the physical point 
can be easily identified as h — (m^/mp) 2 ~ 0.031. Furthermore the entire range of quark mass from to oo can be 
parameterized in a finite interval < h < 1 . Taking advantage of the latter feature, we attempt a global fit including 
the data in the quenched limit f|(||[ . We find that a Pade-type ansatz extending (|l^) given by 

Ppc = cf + cfh^ 1 -^ (26) 
1 + Cp h 

reproduces our data well, as shown in Fig. 15L with x 2 /Nop — 0.33 for Np>p = 5. The fitted value for the chiral 
transition point (3 ct = 1.41(6) is too low, however, presumably because results toward larger values of h including that 
in the quenched limit (h — 1), having small errors, dominantly determine the fit. 



VI. CHIRAL TRANSITION TEMPERATURE 



In order to calculate the transition temperature, we evaluate T pc /my — l/(N t mya) where my is the vector meson 



mass at the transition point (/3, K t (p)) evaluated at zero temperature. The results are plotted in Fig. 16 as a 
function of (mps/rov) 2 - We see that the difference in the values of T pc /my at the chiral transition point (3 c t where 
( m ps/ m v) 2 = and at the physical point (mps/my) 2 = 0.031 is negligible compared with the current magnitude 
of errors. Our values of T pc /my are slightly smaller than those obtained with the standard Wilson quark action on 
iVt = 4 lattices, and consistent with the results from other improved actions Q. 

To estimate the transition temperature in the chiral limit, we interpolate the vector meson mass along the zero- 
temperature critical line K C (T = 0) to (3 — (3 c t. For the numerical value of (3 c t we take the result ( p2|) based on 0(4) 
scaling of the K t line in terms of the Ward-Takahashi identity quark mass. This gives 

T c /m v = 0.2224(51), (27) 

which is plotted by open circle in Fig. |l6]. Converting to physical units, we obtain 

T c = 171(4) MeV (28) 



The solid line drawn in Fig. 16 is a guide to the eyes, representing a global fit of the data to a Pade type ansatz (|2 
with (3 pc replaced by T pc /my and employing the 0(4) values for the exponents. 



VII. CONCLUSIONS 



We have studied the phase structure and the nature of the chiral transition in two-flavor QCD at finite temperatures 
using an RG-improved gauge action and a clover-improved Wilson-type quark action on a lattice with the temporal 
size N t = 4. We have identified the boundary of the parity-broken phase and the finite temperature transition line on 
a two-parameter space of the coupling and the hopping parameter. The chiral transition point is found to be located 
very close to the cusp point of the parity-broken phase. 

A subtracted chiral condensate is shown to satisfy the scaling behavior with the exponents and the scaling function 
universal to the 0(4) Heisenberg model. The quark mass dependence of the transition point is also consistent with the 
0(4) prediction. These results, in agreement with the previous study of the standard Wilson quark action combined 
with the RG-improved gauge action jll]|L7| , indicate that the chiral transition of two-flavor QCD is of second order 
in the continuum limit. 

Assuming the 0(4) scaling, we have extrapolated the transition point towards the chiral limit. Fixing the lattice 
scale in terms of the p meson mass, we obtain T c — 171(4) MeV for the transition temperature. 
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£ 
1.800 
1.825 
1.850 
1.865 
1.875 



K 

0.1300-0.1450 

0.1425 
0.1250-0.1440 

0.1400 
0.1350-0.1400 



traj. 
500-2000 

1000 
500-1900 

2000 
1580-2000 



therm. 
200-500 

300 
200-300 

300 

200 
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1.890 


0.1400 


2200 


300 


1.900 


0.1250-0.1425 


500-2000 


200-400 


1.910 


0.1350 


2100 


300 


1.925 


0.1300-0.1400 


1000-2000 


200-300 


1.950 


n i inn n i A i n 

0.1200-0.1410 


rAA OAAA 

500-2000 


200 


1.975 


0.1300 


1000 


200 


2.000 


0.1150-0.1390 


500-2000 


200-300 


2.100 


0.0900-0.1375 


500-1000 


200-900 


2.200 


0.0700-0.1365 


500 


200 


TABLE I. Simulation parameters on the 16 x 4 lattice. 


a 
P 


A 


traj . 


therm. 


i on 
1.80 


n i onn n i /i ca 


onn 
2U0 


onn rnn 
2UU-OUU 


1 OK 

1.00 


n iocn n i a a n 
U.lzoU— U.144U 


onn onn 
zUO-300 


i nn onn 


i nn 
1.90 


U.lzoU— U.14zo 


onn 
zUO 


onn /inn 


1.95 


n i onn n i a i n 


onn onn 
zUO-300 


i nn a nn 


2.U0 


nil a i onn 


onn onn 
zUO-300 


i nn onn 
1UU-2UU 


o i n 
Z.1U 


a nnnn n i otk 


onn 
300 


onn ititn 


o on 
Z.zU 


a nTnn a 1 o£!c; 


onn onn 
zUO-300 


i nn onn 


2.25 


0.1300-0.1360 


250 






TABLE II. Simulation parameters on the 16 4 (/? = 1. 


8-2.2) and 16 3 x 42 (f3 = 2.25) lattice. 




P 


A 


traj. 


therm. 


1.20 


0.1425-0.1475 


118-120 


70-72 


1.30 


0.1475-0.1525 


118-120 


70-104 


1.40 


0.1500-0.1540 


1 -1 OA 

120-180 


110-280 


1.40 


0.1625-0.1750 


84-100 


36-80 


1.50 


0.1475-0.1540 


100-120 


70 


1.50 


0.1 585-0. 1 700 


88-100 


50-90 


1.60 


0.1400-0.1525 


250-300 


75-510 


1.60 


0.1560-0.1700 


92-500 


50-300 


1.65 


0.1400-0.1500 


300-500 


200-460 


1.65 


0.1 565-0. 1 700 


250-500 


100-400 


1.70 


0.1400-0.1500 


250-300 


150-300 


1.70 


0.1 545-0. 1 700 


250 


80-240 




TABLE III. Simulation parameters 


on the 8 3 x 4 lattice. 




P 


K 


traj. 


therm. 


1.30 


0.1 375-0. 1 500 


120 


55-70 


1.40 


0.1400-0.1525 


90-120 


50-70 


1.50 


0.1375-0.1525 


90-120 


50-80 


1.60 


0.1370-0.1510 


700-1240 


50-150 


1.70 


0.1340-0.1510 


120-1275 


150 



TABLE IV. Simulation parameters on 12 3 x 24 lattice. 
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p 


a' 1 [GeV] 


K C (T = 0) 


Kj. Dwer (iV f = 4) 


Kc PpcT (N t = 


1.20 






0.15044(20) 




1.238(12) 








0.17000 


1.30 


0.563(21) 


0.154082(56) 


0.15475(25) 




1.303(07) 








0.16500 


1.364(04) 








0.16250 


1.40 


0.606(31) 


0.156308(62) 


0.15702(15) 


0.16103(26) 


1.432(05) 








0.16000 


1.50 


0.658(11) 


0.156641(58) 


0.15730(22) 


0.15780(21) 


1.60 


0.707(10) 


0.155259(31) 






1.70 


0.780(10) 


0.151987(22) 






1.80(*) 


0.995(19) 


0.147678(15) 






1.85 


1.033(29) 


0.145526(58) 






1.90 


1.207(36) 


0.143737(48) 






1.95(*) 


1.331(23) 


0.142072(14) 






2.00 


1.445(33) 


0.140811(55) 






2.10(*) 


1.851(58) 


0.139020(21) 






2.20(*) 


2.42(14) 


0.137658(53) 






2.25 


2.40(24) 


0.137225(92) 







TABLE V. Results for the lattice scale a \ the chiral limit K C {T = 0), and the upper and lower boundaries of the 
parity-broken phase K c (N t = 4). Results for K C (T — 0) and a -1 at (5 marked by (*) are obtained using hadron mass results 
from Q. 



p 


K t (N t =4) 


mps/my 


Tpc/mv 


1.600 


0.1543(10) 


0.346(153) 


0.217(11) 


1.650 


0.1533(10) 






1.700 


0.1510(10) 


0.396(170) 


0.234(17) 


1.800 


0.1445(14) 


0.690(92) 


0.211(15) 


1.850 


0.14019(18) 


0.7905(60) 


0.1917(20) 


1.900 


0.13621(15) 


0.8525(39) 


0.1801(12) 


1.925 


0.13417(23) 






1.950 


0.13040(97) 


0.9051(64) 


0.1572(62) 


2.000 


0.12371(73) 


0.9450(36) 


0.1398(29) 


2.100 


0.10921(43) 


0.9790(13) 


0.1114(09) 



TABLE VI. Finite temperature transition/crossover point Kt for N t = 4. Results for mps(T = 0)/mv(T = 0) and 
T pc /my(T = 0) interpolated to the Kt point are also listed. No zero-temperature simulations were made at j3 = 1.65 and 1.925. 
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FIG. 1. Chiral extrapolation of vector meson mass as a function of (mpsa) 2 obtained on 12 3 x 24 lattice. 
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FIG. 2. Phase diagram for the RG-improved gauge action and clover quark action at Nt = 4. 
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FIG. 3. Chiral extrapolation of pseudo scalar meson mass as a function of K on 12 3 x 24 lattice. 
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FIG. 4. Polyakov line obtained on 16 3 x 4 lattice. 
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FIG. 5. Polyakov line susceptibility obtained on 16 3 x 4 lattice. Several data are omitted for clarity of the plot. 
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FIG. 6. Spatial plaquette obtained on 8 3 x 4 lattice. 



14 




1G 



1.5 



1.0 



V 



0.5 



0.0 



0.0 




op=1.80 
□ p=1 .85 
<>P=1.875 
A p=1 .90 
< p=1 .925 
VP=1.95 




0.5 



t/h 



1/p8 



.1.0 



1.5 



FIG. 11. 0(4) scaling fit with h = 2m q a, using all data at f} = 1.8-1.95 and 2m q a < 0.9 on 16 3 x 4 lattice. The best fit is 
obtained at f3 ct = 1.469. 
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FIG. 12. The same as Fig. |ll] using m q a obtained on the zero-temperature 16 4 lattice. The best fit is obtained at f3 c t = 1.462. 
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m q a 

FIG. 13. The pseudo-critical point f3 pc as a function of the Ward identity quark mass m q a. Result of the fit using the 0(4) 
exponent is shown by a solid line. The prediction for (3 c t in the chiral limit from this fit is shown by a star. The location of the 
cusp of the parity-broken phase is shown by a diamond, which is slightly shifted in the horizontal direction for clarity of the 
figure. 
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FIG. 14. The same as Fig. |l| for h = 1/K t - 1/K C . 
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FIG. 16. Chiral transition temperature in units of vector meson mass as a function of (mps/mv) 2 . Hadron masses are 
measured at the same simulation point on a zero-temperature lattice. The chiral transition temperature estimated from (5 c t 
discussed in Sec. M is shown with an open symbol. The solid line is a guide to the eyes based on a Pade-type ansatz. 
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